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Abstract. Chemotactic swimming behaviors of spermato¬ 
zoa toward an egg have been reported in various species. 
The strategies underlying these behaviors, however, are 
poorly understood. We focused on two types of chemotaxis, 
one in the siphonophores and the second in the ascidians, 
and then proposed two models based on experimental data. 
Both models assumed that the radius of the path curvature 
of a swimming spermatozoon depends on [Ca 2 + ],, the in¬ 
tracellular calcium concentration. The chemotaxis in the 
siphonophores could be simulated in a model that assumes 
that [Ca 2 ' ], depends on the local concentration of the at- 
tractant in the vicinity of the spermatozoon and that a 
substantial time period is required for the clearance of 
transient high [Ca 2+ ] t . In the case of ascidians, trajectories 
similar to those in experiments could be adequately simu¬ 
lated by a variant of this model that assumes that [Ca~ + ]j 
depends on the time derivative of the attractant concentra¬ 
tion. The properties of these strategies and future problems 
are discussed in relation to these models. 

Introduction 

After the first discovery of the chemotactic behavior of 
spermatozoa toward an egg by Dan (1950), sperm chemo¬ 
taxis has been reported in various animal species (Miller, 
1977; Eisenbach, 1999), and chemoattractants that mediate 
such behaviors have been identified in a few cases (Ward et 
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al., 1985; Coll et al., 1994; Olson et ah. 2001; Yoshida et 
ah. 2002). Furthermore, recent cell biological studies re¬ 
vealed some of the intracellular enzymatic signaling cas¬ 
cade evoked by the attractant stimulations (Yoshida et ah, 
2003; Kaupp et ah. 2003). However, a major question that 
anyone who observed this phenomenon has to bear in mind 
remains poorly answered: How does such a tiny spermato¬ 
zoon succeed in finding an attractant source? 

What makes this question both more mysterious and 
more interesting is the fact that distinct swimming trajecto¬ 
ries have been observed for spermatozoa of different spe¬ 
cies, which implies that different strategies for chemotaxis 
may underlie these behaviors. For example, the chemotactic 
behavior in the siphonophores has been described as follows 
(Cosson et ah, 1984): spermatozoa show trajectories of 
large diameter (700-1000 jam) while swimming far from 
the “cupule," which in these species is an extracellular 
structure of the egg and serves as an attractant source, and 
trajectories of smaller diameter (200 jam) in the vicinity of 
the cupule; the transition between the two modes is progres¬ 
sive. In the ascidians, on the other hand, spermatozoa ex¬ 
hibit characteristic “chemotactic turns’' during chemotaxis, 
but the curvature of the trajectories has no noticeable de¬ 
pendence on the distance to the attractant source (Yoshida et 
ah, 1993, 2002; lshikawa, 2000). The behavior of sea urchin 
spermatozoa seems to be a little more complicated but still 
exhibits unambiguous chemotaxis (Ward et ah, 1985). 
Thus, it appears that spermatozoa in different species often 
respond to the attractant differently, although the goal is the 
same in all cases—to find an egg and facilitate fertilization. 

One approach to revealing the underlying strategies for 
sperm chemotaxis may be to build a model based on the 
experimental observations and then study how the model 
works. In the present paper, we propose two models, one 
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applied to the siphonophores and the second to the ascid- 
ians. We focused on these two taxa because their swimming 
behaviors are simple but clearly distinct, and also because 
substantial experimental data are available in the literature 
(see Results for details). These experimental results were 
simplified, as summarized in Figure I, and incorporated into 
the models, which required only a few simple and reason¬ 
able assumptions. We show that numerical calculations of 
sperm trajectory using the two models result in trajectories 
very similar to those observed experimentally. 

Materials and Methods 

General remarks on the models 

It has been suggested that spermatozoa swim with three- 
dimensional beating waves whose characteristics coniine 
them to a two-dimensional space when they are confronted 
with the interfaces between two media such as water/air or 
glass/water (Cosson et ai, 2003). For this reason as well as 
to simplify the observation of sperm trajectories under a 
microscope, our experiments have been restricted mainly to 


spermatozoa located at the glass/water interface of the lab¬ 
oratory chambers. Since experimental background is neces¬ 
sary to test the validity of a model, the models we used were 
also limited to two dimensions. At the interface, most sper¬ 
matozoa of siphonophores and ascidians show circular 
movement, usually in one preferred direction. A reasonable 
explanation is that the movement in the three-dimensional 
free environment is helical and that some of the initially 
homogeneously suspended spermatozoa were going upward 
and tended to be confined to the air/water interface, while 
others were going downward, to the water/glass interface, as 
was shown for sea urchin spermatozoa (see Cosson et ai, 
2003). Sperm confined to the glass surface still exhibit 
unambiguous chemotaxis as reported thus far, and this be¬ 
havior is our topic in the present study. 

In our model, a spermatozoon is treated as a mathematical 
point. This point circulates in the .v-y plane in one direction 
(clockwise for the ascidian spermatozoa observed in the 
vicinity of the glass slide surface, and counterclockwise for 
the siphonophore spermatozoa), with a radius of curvature 
that changes depending on the attractant concentration and 
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Figure 1. A summary, based on previous reports, of spermatozoa behavior in the presence of the attractant 
in two species of siphonophore (Muggiaea kochi, Chelopheyes appnulicuiata) and one species of ascidian 
i'iona iiitesthuilis). (A) tn the siphonophores. the curvature is negatively dependent on the attractant concen¬ 
tration in artificial seawater (ASW. top), hut the dependence is lost in Ca 2+ -free seawater (CFSW, bottom). (B) 
In the ascidians. no such dependence is observed, even in ASW (top), but a rapid modulation of the curvature 
leading 1 a turning behavior occurs in the presence of the attractant gradient (bottom left). Arrows indicate 
direction oi the attractant source. The curvature is larger when swimming toward the source and smaller when 
swimming n the opposite direction. A short delay in the response has been often observed in this correlation 
(Ishikawa. 2000: Yoshida et oi. 2002). In CFSW, no such modulation in the curvature is observed, even under 
the attractant gradient (bottom left). 
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on time. Let us define some parameters: t , time; pit), 
position of sperm at /; p(f) = v{cos (c pit)), sin i(f)it))}. 
velocity vector of sperm; rit), curvature; and dp), attract- 
ant concentration sensed by the sperm at position p (Fig. 
2A). Since no significant changes in the swimming velocity 
were found during the chemotactic behavior in these species 
(Cosson et ah, 1984; Ishikawa, 2000; Yoshida et ah , 2002), 
let us assume the amplitude of the velocity vector as con¬ 
stant ( = v). Numerical integration of sperm movement was 
evaluated as follows with a time step (A/) of 1.0 ms: 

pit + Ar) = pit) + AM/) (Eq. 1) 

4>(t + A/) = 4>it) + Atv/rit) (Eq. 2) 

Updating the angular component of velocity, but not the 
vector itself, as in the equation (2), moderates the accumu- 
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Figure 2. (A) A schematic drawing that shows the four parameters 

used in the calculation of sperm trajectory. O, the origin; p, position at t 
(vector); v, velocity (vector); (/>, angular component of v (scalar); r, radius 
of curvature (scalar). (B) A computed snapshot of the chemoattractant 
gradient used in the present study. The attractant concentration profile in 
the square area of ± 1 mm of the x-y plane is plotted in the c-axis. The peak 
(hilltop) corresponds to the chemoattractant source located at the origin 
(x - y = 0). 


lation of error during the integrations for a long time period. 
All the calculations and graphical outputs were done by 
using the Mathematica 3.0 program (Wolfram Research). 

Attractant profile 

For most observations of chemotaxis, the following ex¬ 
perimental design was adopted: a glass pipette filled with 
chemoattractant or any fluid (within agar-gel in many cases) 
to be tested was placed in a drop of water containing 
dispersed sperm cells. This allows a spatial gradient of 
attractant concentration to be established rapidly: sperm 
located in the vicinity of an attractant source usually ap¬ 
proach it within no more than tens of seconds. Since mol¬ 
ecules such as an organic compound or a small protein 
normally have diffusion coefficients in water of ~10 10 
nr/s or smaller, it is an acceptable approximation that the 
attractant profile does not change much during the approach 
of the spermatozoon. In our model, therefore, we simply 
used a computed snapshot of a solution of a diffusion 
equation as an attractant profile (Fig. 2B). Moreover, the 
drop is largely spread on the slide but is much thinner in the 
^-axis, so variation of attractant concentration along the 
:-axis should be small. Therefore, we consider the attractant 
gradient field with a two-dimensional diffusion equation 
with a coefficient of D: 

8c/8t = D(8 2 c/8.\ 2 + 8 2 c/8y 2 ) (Eq. 3) 

In the polar coordinates of p, 0 (.v = p cos (0), v = 
p sin (0)), equation (3) can be expressed as 

8c/8t = D(8 2 c/8p 2 + l/p5c/5p + l/p 2 S 2 c/S0 2 ) (Eq. 4) 

Since we assume a situation of spatial symmetry, the last 
term of equation (4) can be regarded to be zero. 

8c/8t = Di 8 2 c/8p 2 + l/p5c/Sp) (Eq. 5) 

As a boundary condition, we assumed that the concentration 
at the origin (/.<?., the tip of a pipette) is always kept at a high 
level (c(0) = 1). We set D as 2.0 X 10 10 nrr/s and used 
a numerical solution of equation (5) at / = 60 s as the 
profile for subsequent trajectory calculations (Fig. 2B). We 
also tried some other static attractant profile models (expo¬ 
nential, parabolic, solution of the diffusion equation at dif¬ 
ferent times, etc.) and found that most definitions resulted in 
qualitatively similar responses to what we show next. 

Results 


The siphonophore model 

Some properties of sperm-swimming in siphonophores 
are summarized in Figure 1A. This summary is based on 
observations of Muggiaea kochi and Chelopheyes appen- 
dicukita (Cosson et ah, 1984). In artificial seawater (ASW), 
the spermatozoa show trajectories of larger diameter (700- 
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Figure 3. Simulations of the siphonophorc model. (A) A cubic tunc 
tion as one of the examples ol a function that negatively correlates 
attractanl concentration to curvature c?t the trajectory. (B) A trajectory 


1000 /xm) while they are far from the attractant source (—5 
mm), and smaller trajectories (—200 /xm) near the attractant 
source (—0.2 mm). However, this dependence on distance 
from the attractant was lost in the absence of external Ca 2 ", 
which suggests that Ca 2+ influx is involved in the modula¬ 
tion of the curvature of trajectories. We incorporated these 
findings into the assumptions set in the model as follows: 

(a) Curvature decreases as intracellular Ca 2 ' concentra¬ 
tion, ([Ca 2+ ],), increases ( i.e ., a negative correlation). 

(b) [Ca 2 + ], positively correlates with c(p), the attractant 
concentration which is sensed by the sperm. 

Assumption (a) also agrees with the evidence that sper¬ 
matozoa treated with a Ca 2 ^ ionophore showed a curvature 
of —200 /xni, which increased up to —800 /xm at lower 
Ca 2+ concentration. Since any quantitative relationship 
among these factors in real spermatozoa is not yet available, 
we set some trial functions and studied the behavior. We 
first tried a simple linear correlation between [Ca 2+ ], and 
dp), and a cubic function for [Ca 2 ']j and r(t) (Fig. 3A): 

[Ca :+ ], = c(p) (0 < c < 1; [arbitral unit (a.u.)]) (Eq. 6) 

i (/) = 400(1 - [Ca 2 l? [/xm] (0 < [Ca 2 ] < 1; [a.u.]) 

(Eq. 7) 

This now allows us to study how spermatozoa swim under 
these conditions. The results show that the trajectories gen¬ 
erated by such conditions are far from cheniotactic behav¬ 
ior. Spermatozoa do circle around the attractant source but 
never finally approach the egg (Fig. 3B). This result is not 
surprising because our present system is fully symmetric in 
terms of time and space components. When using functions 
other than those of equations (6 and 7), no chemotaxis is 
observed at all. However, a remarkable change occurs if we 
incorporate an additional condition: 

(c) [Ca 2 + 1 ; has different rates of increase and decrease 
(slow to decrease). 

This condition is not unlikely, because pumping of the 
cytosolic Ca 2+ out of the cell or into the intracellular store 
against the electrochemical gradient across the membrane is 
an enzymatic process, which is energy-consuming, whereas 
Ca 2+ influx generally occurs rapidly. Therefore, we then 
incorporated this additional condition as follows: 

A[Ca 2 + ],/\t = - l/T([Ca 2+ ], - dp)) 

when [Ca 2+ ],>dp) (Eq. 6') 


simulated with the function shown in (A). Time to lower (Ca 2 "], is not 
considered. The sperm started at position (.v. y) = ( -800 jam, -800 /am) 
with an initial angle of 0.1 radian and a constant speed ot 500 /am/s 
(arrow). The trajectory during the initial 40 s is shown. (C) A trajectory 
simulated under the same conditions as in (B), except that t ot 0.5 s is 
incorporated as a decay-time constant to lower [Ca 2+ ],- 
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Figure 4. Simulations of the ascidian model. (A) A sigmoid function 
as one of the examples of a function that positively correlates time-derived 


where r is the decay time constant. This condition was 
applied only when [Ca 2 |, > c(p)\ otherwise, equation (6) 
was applied as before. When trajectories were calculated 
with these new sets of conditions (6, 6', 7), we found that 
the incorporation of a time component for the Ca 2+ decay 
results in chemotactic trajectories (Fig. 3C) that are not 
unlike those obtained in the experimental observations. 
Likewise, chemotactic trajectories were reconstructed for 
a wide range of arbitrarily selected functions that relate 
"dp)' to ‘[Ca - ],* positively instead of by equation (6), 
and relate ‘[Ca 2 ‘ [,’ to h r(t)' negatively instead of by 
equation (7), even though some of them, of course, showed 
heavily distorted trajectories of approach toward the 
attractant source and unusually long or short times of 
approach. 


The ciscidhm model 

The behavior of ascidian sperm is quite different from 
that in the siphonophore sperm. The behavior illustrated 
in Figure IB is based on observations of Ciona intesti¬ 
nal is by lshikawa (2000) and Yoshida et al . (1993). In the 
presence of an attractant without a gradient, the curvature 
of the sperm trajectory is independent of the concentra¬ 
tion of the attractant. and is close to that in the absence 
of the attractant (Fig. IB). Under the gradient of attract¬ 
ant concentration, however, spermatozoa show rapid 
changes in their track diameters, so-called chemotactic 
turns. Since this turning behavior is lost in the calcium- 
free seawater (CFSW) or in the presence of Ca 2 ~ channel 
inhibitors, it has been suggested that rapid changes of 
diameter are dependent on [Ca 2 + ] i . Therefore, we again 
incorporated condition (a), as we did in the siphonophore 
model. Next, we need to infer the relationship between 
[Ca 2 ]j and c[p), the attractant concentration that the 
spermatozoan senses. Since modulation of diameter oc¬ 
curs only in the presence of an attractant gradient, we 
assumed that [Ca 24- ], depends on the temporal changes of 
tip) but not on the “absolute" concentration itself. Thus, 
we assume a second condition: 

(d) [Ca - " ^ negatively correlates to the time derivative of 
the attractant concentration. 

With conditions (a) and (d), it follows that r(t) positively 
correlates to Ac{p(t))/At, the time derivative of the attract- 


attractant concentration to the curvature of the trajectory. (B) A trajectory 
simulated with the function shown in (A). The sperm started at position (.v, 
y) = (-750 ix m. 750 jam) with an initial angle of 0.8 radian and a 
constant speed of 250 jum/s (arrow). The trajectory during the initial 40 s 
is shown. (C) A trajectory simulated under the same conditions as in (B), 
except that a delay of 150 ms in the response of curvature to the time- 
derived attractant concentration is incorporated. 
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ant concentration. This assumption was incorporated in the 
ascidian model as follows: 

r(t) = F(±c(p(t))/\t) (Eq. 8) 

Where F is the function that positively correlates be/St with 
r(t). We studied sperm behavior with a sigmoid function for 
F (Fig. 4A) and found that this condition alone is sufficient 
to show chemotactic trajectories with successive turns when 
approaching the attractant source (Fig. 4B). We further 
incorporated a delay in the sperm response of r(t) to 
\c(p(t))/St with a time range of tens of milliseconds, 
because such a delay has been found in the analysis of the 
trajectory of real sperm (Ishikawa, 2000; Yoshida et al, 
2002). This resulted in similar trajectories, but with a twist 
(Fig. 4C), which is often observed in real chemotaxis in the 
ascidian. Thus, we find that the delay of sperm response is 
not an absolute necessity for chemotactic behavior in the 
case of the ascidian model, but this parameter results in 
trajectories that are somewhat more realistic. 

Discussion 

In the present study, we proposed two comprehensive 
models for strategies of sperm chemotaxis in the siphono- 
phores and the ascidians. With these models, chemotactic 
trajectories similar to those observed for real spermatozoa 
could be reconstructed. We found that there are at least two 
ways to identify the location of the “hilltop"—the source of 
the chemoattractant—without looking around for it. The 
siphonophore's way is to “walk" circularly, with large- 
diameter curvatures at low altitude and smaller ones at 
higher altitude. This alone is not enough for success; how¬ 
ever, success can be achieved if any increase in curvature is 
followed by a time lag. In the ascidian’s way, one needs to 
sense the steepness (gradient) but not the height (absolute 
concentration): the curvature is large when climbing steeply 
uphill, medium when the slope is gentle, and small when 


going downhill. The time-delay condition is not always 
necessary, but such a delay introduces twists into the tra¬ 
jectory. 

Even though these models employ some conditions based 
on experimental results, these conditions may seem to be 
oversimplified. However, the goal of the present modeling 
study is not to simulate the natural behavior of the sperma¬ 
tozoa perfectly but to find out what elements are needed to 
reconstruct the phenomenon of interest. For this purpose, it 
is necessary to focus on a small number of important pa¬ 
rameters. 

Nonetheless, some lines of evidence in addition to those 
mentioned earlier, support our simplification of the condi¬ 
tions. First, it has been shown that asymmetrical bending 
waves are induced when a high concentration of Ca 2 ~ is 
applied to the flagella of demembranated sea urchin sperm 
(Brokaw, 1979). This may support our condition (a). Next, 
it was recently suggested that a store-operated Ca 2+ channel 
(SOC) regulates chemotaxis in the ascidian (Yoshida et al, 
2002). Involvement of such a SOC may account for the time 
delay in the curvature response to the temporal changes in 
the attractant concentration: this feature was often observed 
in the experiments and is introduced into the ascidian model 
(Fig. 4C), since activation of the SOC requires depletion of 
internal Ca 2+ stores evoked by the releasing signals such as 
inositol 1,4,5-phosphate. 

Quantitative comparison of the model outputs with ex¬ 
perimental results is very important for validation of the 
models. Unfortunately, we are not ready to do such a study 
because the models currently have high degrees of freedom. 
The two functions, one that relates attractant stimulus with 
[Ca 2 + ], and the other that relates [Ca 2+ ], with curvature, 
cannot yet be quantitatively defined. Of course, one desir¬ 
able experiment might be to experimentally define these 
functions and then to carry out a quantitative comparison. 
But this requires many very difficult technical break- 



Figure 5. Changes in [Ca 2+ ], during chemotaxis. (A) Plot of [Ca 2+ ], as a function of time in the modeled 
chemotaxis for siphonophores (Fig. 3C). (B) Plot of [Ca 2+ ], versus lime in the modeled chemotaxis tor the 
ascidian d ig. 4B). It is assumed that [Ca 2+ |, is proportional to the inverse of the radius of curvature. 
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Figure 6 . A difference in the sensitivity to a local attractant peak 
between the siphonophore and ascidian strategies. (A) A secondary source 
with one-fourth the attractant concentration of the main one is added at the 
position (.v. v) = (-500 /am, -500 /am). (B) Sperm trajectory with the 
siphonophore model under the dual peak profile in (A). The trajectory is 
superimposed upon the attractant profile. Note that the spermatozoa found 
the main peak. The initial condition is same as that in Fig. 3B. (C) Sperm 
trajectory simulated with the ascidian model. Note that the spermatozoa 
now found the local peak instead of the main peak. The initial condition is 
the same as that in Fig. 4B. Arrows in (B) and (C) indicate the initial 
positions and directions. 

throughs such as quantification of local attractant concen¬ 
tration, application of an attractant stimulus varying with 
time, and so on. Another indirect but more practical exper¬ 


iment would be to measure [Ca 2+ ], during chemotaxis. Our 
models predict that there will be a temporal [Ca 2+ ] 1 pattern 
during chemotaxis specific to each model: [Ca 2+ ], oscillates 
in both models, and the base level elevates in the siphono- 
phores but not in the ascidians (Fig. 5). Even though mea¬ 
surement of [Ca 24 ]; in swimming spermatozoa is still tech¬ 
nically challenging, mainly due to the small volume of 
cytoplasm and fast movements of the sperm, we hope that 
recent innovations in Ca 2 indicators and image sensors 
will make it possible in the future. 

What is the significance of the difference in the two 
strategies? Since temporal change of the attractant concen¬ 
tration is the critical parameter for the modulation of cur¬ 
vature in the ascidian model, we expect that the ascidian 
spermatozoon is more sensitive to local change of concen¬ 
tration than the siphonophore spermatozoon, in which the 
absolute concentration value is the most critical. Simulation 
under two attractant sources, one of high concentration and 
the other of low concentration, led to the expected results: 
the ascidian spermatozoa found the lower peak of attractant 
in the vicinity, and the siphonophore spermatozoa reached 
the higher peak (Fig. 6). Thus, an interesting future problem 
would be to test this prediction in biological experiments. 
One may also be interested in possible relationships be¬ 
tween such characteristics of sperm strategy and the envi¬ 
ronment that the species are facing. However, such a ques¬ 
tion cannot be addressed currently. Most of our knowledge 
has been limited to sperm in two-dimensional laboratory 
conditions, even though the goal should be to reach a deep 
understanding of strategy in the three-dimensional natural 
environment. Since theoretical treatments suggests that drag 
forces near the interface have substantial effects on the 
flagellar motion of sperm (Katz, 1974), we should avoid 
applying knowledge gained in two dimensions thought¬ 
lessly in building three-dimensional models. One needs to 
know how parameters that describe three-dimensional heli¬ 
cal motions are affected by the attractant gradient. To ac¬ 
complish this goal, we hope that techniques to measure 
sperm trajectories in the three-dimensional environment, as 
pioneered by Crenshaw et al. (2000). will soon be more 
accessible. 
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